Elastic scattering losses in the four-wave mixing of Bose Einstein Condensates 
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We introduce a classical stochastic field method that accounts for the quantum fluctuations re- 
sponsible for spontaneous initiation of various atom optics processes. We assume a delta-correlated 
Gaussian noise in all initially empty modes of atomic field. Its strength is determined by comparison 
with the analytical results for two colliding condensates in the low loss limit. Our method is applied 
to the atomic four wave mixing experiment performed at MIT [Vogels et. al, Phys. Rev. Lett. 89, 
020401, (2002)], for the first time reproducing experimental data. 
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In recent years we observe a growing number of ex- 
periments in which the atomic Bose-Einstein condensate 
evolves in a nontrivial way. A whole new area of non- 
linear atom optics was born. The most striking example 
of such a nonlinear process is the atomic four- wave mix- 
ing (4WM). In close analogy to its optical counterpart, 
atomic 4WM consists of generation of the new atomic 
beam in the nonlinear interaction of three overlapping 
matter waves. For the main part atomic 4WM is an ex- 
ample of a stimulated process. However, during this pro- 
cess there are also collisions between individual atoms 
that lead to a population of initially unoccupied atomic 
states. These processes have spontaneous initiation but 
by nature they are also examples of the four particle mix- 
ing. They amount to the elastic scattering losses from the 
coherently evolving condensates. 

The standard tool used to describe dynamics of the 
condensate within mean field approximation is the cele- 
brated Gross-Pitaevskii equation (GPE). As it stands, it 
is capable of describing stimulated processes but not the 
spontaneous ones. However, at least in some experiments 
[l| the elastic scattering losses may become significant. 
There are at least two theoretical attempts to estimate 
such losses during the collision of condensates. In the first 
one the authors used momentum-dependent higher 
order correction to the nonlinear coupling constant in 
the GPE, introducing complex scattering length. In the 
second one ^ Q [1| the field theoretical formulation was 
used. To make it effective, the authors approximate the 
second quantized hamiltonian by a quadratic form. Both 
methods give very similar results but are applicable only 
if the elastic scattering losses are merely a small correc- 
tion. 

It is the purpose of this Letter to formulate a gen- 
eral method of describing elastic scattering losses in the 
nonlinear atom optics processes. To this end we add to 
the GPE a classical gaussian noise, representing vacuum 
quantum fluctuations of the atomic field and an auxil- 
iary field holding atoms scattered out from BEG. Such a 
technique has its roots in quantum optics. 

Spontaneous optical processes have their origins in 



quantum fluctuations. The best known example is a 
process of superfluorescence 0. In this phenomenon a 
sample of atoms is prepared in the internal excited state. 
Spontaneously emitted photons create an avalanche of 
photon emissions. When the light field becomes strong, 
it is well described by a classical electromagnetic field. 
However the initiation has a quantal nature. This quan- 
tum initiation was successfully imitated by a classical 
noise 7:]. There are also general methods of map- 
ping quantum fluctuations into stochastic term in the 
evolution equations of quantum optics (generalised P- 
representation methods) 0. 

In optics one can find numerous other processes initi- 
ated by spontaneous emission and eventually upon pop- 
ulating empty modes turning into stimulated processes; 
eg. spontaneous Raman scattering 9], parametric down 
conversion .IQj . etc. There are also similar examples in 
atomic and molecular physics llj. Our method is general 
and is capable of treating many of these processes. Here 
we demonstrate the method using the 4WM of coherent 
matter waves. 

The first experiment demonstrating 4WM in a sodium 
Bose-Einstein condensate was performed at NIST 12j. 
This was followed by a theoretical and numerical treat- 
ment of the experiment |l3L Il4l|. In the experiment, a 
short time of free expansion of the condensate, after it 
was released from the magnetic trap, was followed by 
a set of two Bragg pulses 0|, which created moving 
wavepackets. These wavepackets, together with the re- 
maining stationary condensate, due to nonlinear interac- 
tion and under phase matching conditions created a new 
momentum component in the 4WM process. The stan- 
dard starting point for the description of atomic 4WM 
process is the Gross-Pitaevskii equation 
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Here is the total number of atoms, l^tp is proportional 
to the atomic number density and is normalized to one, 
g = 'inh^a/m is the nonlinear interaction strength, m 
the atomic mass, a is the scattering length and is a 
confining potential. A compact ground state wavefunc- 



tion 4'(f, 0) is created in harmonic trap potential V and 
centered around r = with \E'(0, 0) = the maximum 
value. Once this ground state is created, V is turned off. 
The development of ^{f,t) is now described by Eq. 
with V — 0. Later on, a set of Bragg pulses is applied and 
parts of the condensate begin to move. We can define two 
timescales characterizing evolution of the condensate: a 
nonlinear interaction time tnl — {gN\'^m\'^ /h)^^ and a 
collision time TcoI- The latter is defined as a time it takes 
two wavepackets uniformly moving along x to move apart 
(so they just touch and cease to overlap), Tcoi = 2rTp/v 
where ttf is the initial radius of the condensate in the 
X direction (Thomas - Fermi approximation), and v is 
the relative velocity. The ratio of these two timescales 
determines the output of the 4WM process. 

The initial condition immediately after application of 
the Bragg pulses at ti, can be approximated as being a 
composition of the BEG wavepackets, identical in shape 
to ^'(r, ii) (for more details see 0|): 

$(f,ti) = vl/(f,ii)^/VV^.--A. (2) 

Here fi — Ni/N is the fraction of atoms in the i-th 
wavepacket and X]i=i fi = ^- ^ ^^^w wavepacket with 
P4 = Pi — P2 + P3 will build up, thanks to the non- 
linear interactions accounted for by the last term in the 
Gross Pitaevskii equation After a while, the fourth 
wavepacket will grow to the macroscopic level. Using the 
de Broglie relations: ki — Pi/Ti and oji — hkf/2m we 
have: 
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= ^$,(f,t)e»(^^--~-.*), (3) 

1=1 

with initial conditions 

Mr, h) = ///'*(r, h), z = 1, 2, 3; ^^(f, t^) = 0. 

Variation of the <i>i is assumed to be slow as compared to 
that of the exponential in equation (|3J. Four equations 
for this slow dependence are obtained from and 
We have, when V is turned off 0: 

Zm ^ — ' 

+ 257VV,+i^*+2V'«+3 (4) 

where we use the convention in which all indices are taken 
modulo 4. To account for the elastically scattered atoms 
we introduce an additional component of the wavefunc- 
tion ipB- It is this part of the wavefunction which will be 
initiated by the classical stochastic field. The stochastic 
field (r, t) must be added to the equation of motion 
Q to trigger the elastic scattering process of two parti- 
cles from the condensate wavefunctions V'i and ijjj to the 
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FIG. 1: Number of elastically scattered atoms from the pair 
of counter-propagating condensate Gaussian wavefunctions, 
with relative velocity 1.75 mm/ s, as a function of time. The 
lower curve corresponds to small number of scattered atoms 
{N = 10^ and gaussian width a = 9.1 nm), when the bosonic 
stimulation does not occur. In this limit all three methods: 
complex scattering length calculation , field theory , and 
our stochastic method give indistinguishable results. The up- 
per pair of curves corresponds to higher number of atoms 
(A^ = 1.6 ■ 10^ and gaussian width a — 15.8 fim) - dashed line 
was obtained using complex scattering length method and the 
solid line is a solution of eq.® 

background field ^b- It is a four particle process and it 
must be implemented for each pair of colliding wavefunc- 
tions i,j in such a way that the total number of atoms 
in colliding waves -I- background field is still conserved. 
The resulting set of equations reads 

{thdt + = <77V(|^,|2 + 2V|V,f + 2|VBnV'. 

2m 

+ 257VV',+iV'.+3^r+2 (5) 
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For numerical calculations we assume that (r , t) is a 
gaussian stochastic process with zero mean and the only 
nonvanishing second order correlation function equal to 
{£.ij{r,t)£_ij{'r' ,t')) = Aij5p^p'6t,t' ■ Here Kronecker delta 
functions are assumed both in space and time since we 
refer to numerical simulations with spacial grid and dis- 
crete time steps. Notice that we assign different stochas- 
tic process to each pair of colliding wavepackets, antici- 
pating dependence on parameters like relative velocity. 

Equations H5I6|) may be obtained from multiatom sys- 
tem hamiltonian upon using Bogoliubov decomposition 



of the atomic field operators into condensate parts and 
initially empty modes ip = XliLi ''Pi + ''Pb in a way anal- 
ogous to that presented in jjj- We do not however ex- 
plicitly decompose -tps into plane waves but obtain it's 
Heisenberg equation of motion assuming only that it 
commutes with operators of macroscopically occupied 
modes i/'i- Finally stochastic field is added to tpB in the 
source terms for 4WM as shown above. These stochas- 
tic terms mimic vacuum quantum fluctuations leading 
to spontaneous elastic scattering loss. But as elastically 
scattered atoms reside inipB they may eventually get am- 
plified via bosonic stimulation when population becomes 
significant. This is an outline of our stochastic method; 
details will be presented elsewhere [f^ . 

To fully determine equations H5I6|I we need to specify 
the value of constants Aij . Just as it has been done in the 
case of superfluorescence mentioned above, we can find 
Aij by the requirement that it reproduces known limiting 
analytic results . In reference the elastic scattering 
losses were computed analytically in perturbative regime 
for two colliding gaussian-shaped wave packets. Further- 
more these wavepackets were assumed to evolve without 
losses and without spreading. The number of elastically 
scattered atoms as a function of time was found in the 
form 

^<"^(^)"-(^'). 

where tr is a width of the gaussian wave-packets and Q 
is the wave vector corresponding to the absolute value of 
the momentum of each of the wave-packets in the cen- 
ter of mass frame. The same quantity might be calcu- 
lated approximately under analogous assumptions using 
the stochastic classical noise. The equation for ipB in this 
case reads 

ihdti^B = -— V2^B+5A^(|V'B|' + 2|Vl|' + 2|V2nV'B 

+ 2gN^iM^B (8) 

where ipi,2 are two counter-propagating gaussian wave- 
functions. The approximations of [3 amounts in retain- 
ing on the right hand side of equation |S1 only the last 
term. The approximate solution obtained this way gives 
the number of elastically scattered atoms as a function 
of time in the form 

Sstochit) = AuAt—— — Erf ^-^t 9 
m[2Q) \ a J \ ma I 

Comparing with ® we obtain Au = j^^- As 
we anticipated A12 depends on the relative velocity of 
wavepackets ipi and ip2, which in our case is equal (2Q). 
Note that once Aij 's are determined our numerical ap- 
proach has no more adjustable parameters. In Fig. ^ 
we are comparing the solution of Q with a numerical 
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FIG. 2: Population of the third wavepacket normalized to the 
initial seed population as a function of time (continuous line) . 
Parameters used in the simulation correspond to experiment 
performed in the Ketterle group Q. Total number of atoms 
equal to 5mln. Also shown: circles - experimentally measured 
values, and dashed and dot lines - solutions of the Gross- 
Pitaevskii equation with real and complex scattering length 
respectively. 



solution of full equation ijSJi. Note the growing discrep- 
ancy between perturbative and non-perturbative results 
for larger losses. They result from bosonic enhancement 
present in the non-perturbative regime. We also stress 
that in the non-perturbative regime the stochastic noise 
is crucial at the early stage of evolution. It may even 
be dropped from equation ((SJ when Bosc enhancement 
takes effect. This is why the strength of classical noise 
may be determined in the perturbative regime. Finally, 
we point out that some analogies regarding the break 
down of perturbative approach were found in the study 
of atom-molecule conversion within positive-P represen- 
tation [rzlliit. 

With the equations H5I6(I fully determined we turn our 
attention to the recent experiment from MIT Q]. This 
experiment, due to the large value of the ratio of collision 
to nonlinear timescales, had very large number of elasti- 
cally scattered atoms. Experimental configuration con- 
sists of two initial wavepackets of equal strength (« N j2 
atoms in each) and the third wavepacket of just a tiny 
fraction of N. Magnetic trap used to generate the Sodium 
condensate had frequencies of 80, 80 and 20 Hz in axial 
direction, hence it has a shape of a cigar. Applied optical 
Bragg pulses to create moving wavepackets propagated 
approximately at the same angle of ~ 0.5 rad with re- 
spect to the long axis of the condensate corresponding 
to a relative velocity of 20 mm/s. In two series of 4WM 
measurements chemical potential of the condensate was 
2.2 and 4.4 kHz, which we identified as corresponding 
to 5 and 30 million atoms respectively. In Figure El we 
plot the population of the third wave-packet (the one 
that was initially seeded) as a function of time. The cir- 
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FIG. 3: Population of the third wavepacket normahzed to 
the initial seed population as a function of time (continuous 
line). Parameters used in the simulation correspond to exper- 
iment performed in the Ketterle group jj. Total number of 
atoms equal to 30 mln. Also shown: circles - experimentally 
measured values, and dashed and dot lines - solutions of the 
Gross-Pitaevskii equation with real and complex scattering 
length respectively. 

cles are the experimental data extracted from paper Q . 
Several theoretical curves are plotted. The dashed line 
represents results neglecting all elastic scattering losses. 
The dotted line accounts for the losses by means of the 
complex scattering length . We see that neither of the 
curves reproduces experimental results. Our stochastic 
method gives the solid line which is much closer to the 
experimental data. It has been computed with the pa- 
rameters of the experiment including the initial number 
of atoms infered from the paper as being equal to 5mln. 

In FigEl similar comparison is made for larger sam- 
ple of 30mln atoms. Again our results reproduce the 
experimental data very well. We feel that remaining dis- 
crepancy (our results seem to be consistently under ex- 
perimental points) is due to indistinguishability of BEG 
and thermal atoms in the region of the momentum space 
occupied by BEG. 

In conclusion: We have formulated the classical 
stochastic field method that accounts for the quantum 
fluctuations responsible for spontaneous initiation of var- 
ious atom optics processes. For instance we can treat 
oscillations between atomic and molecular condensates 
triggered by optical or magnetic field effects 0,0. The 
method is then applied to the atomic 4WM. It gives for 
the first time excellent agreement with the recent MIT 
experiment, where the scattering losses where substan- 
tial. 
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